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Abstract 



■ We present an algorithm to compute the annihilator of (i.e., the linear differential 

equations for) the multi- valued analytic function / A (log f) m in the Weyl algebra D n 
for a given non-constant polynomial /, a non-negative integer m, and a complex 
number A. This algorithm essentially consists in the differentiation with respect to s 
of the annihilator of f s in the ring -D n [s] and ideal quotient computation in D n . The 
obtained differential equations constitute what is called a holonomic system in D- 
module theory. Hence combined with the integration algorithm for D-modules, this 
ON ■ enables us to compute a holonomic system for the integral of a function involving 

, the logarithm of a polynomial with respect to some variables. 

1 Introduction 

For a given function u, it is an interesting problem both in theory and in practice to 
determine the differential equations which u satisfies. Let us restrict our attention to linear 
differential equations with polynomial coefficients. Then our problem can be formulated as 
follows: Let D n be the Weyl algebra, i.e., the ring of differential operators with polynomial 
coefficients in the variables x = (xx, . . . , x n ). An element P of D n is expressed as a finite 
sum 

P= a ^ x ° d ^ (!) 

with x a = i" 1 ---!"", d 13 = d^-"-d^ n and a a> p G C, where a = (a±, . . . ,a n ), (3 = 
(fix, . . . , fi n ) G ~N n are multi-indices with N = {0, 1, 2, ... } and <9; = d/dxi (i — 1, . . . ,n) 
denote derivations. The annihilator of u (in D n ) is denned to be 

Ann DnM = {Pe D n | Pu = 0}, 

which is a left ideal of D n . Since D n is a non-commutative Noetherian ring, there exist 
a finite number of operators Pi, ... , P/v G D n which generate Ann^w as left ideal. Thus 
we can regard the system 

p lU = ■■■ = P N u = 
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of linear (partial or ordinary) differential equations as a maximal one that u satisfies. 

As to systems of linear differential equations, there is a notion of holonomicity, or being 
holonomic, which plays a central role in .D-module theory. See Appendix for a precise 
definition. A holonomic system of linear differential equations admits only a finite number 
of linearly independent solutions although it is not a sufficient condition for holonomicity. 
A holonomic function is by definition a function which satisfies a holonomic system. 

The importance of the holonomicity lies in, in addition to the finiteness property 
above, the fact the it is preserved under basic operations on functions such as sum, 
product, restriction and integration. Hence starting from some basic holonomic functions 
we can construct various holonomic functions by using such operations. 

As one of basic holonomic functions, let us consider f x with a non-constant polynomial 
/ in x — (xi, . . . ,x n ) and a complex number A. Then the function f x is holonomic and 
there is an algorithm to compute its annihilator strictly ([9].[3].[T3]). 

Our purpose is to give an algorithm to compute the annihilator of / A (log/) m with 
a positive integer m and to prove that it is a holonomic function. This is achieved by 
differentiation with respect to the parameter s of the annihilator of f s in D n [s]. This 
method can be extended to functions of the form f Xl ■ ■ ■ /^(log fi) mi ■ ■ ■ (log fN) mN for 
polynomials fk, complex numbers and nonnegative integers m^. 

Since the algorithm yields a holonomic system, we can apply the integration algorithm 
for .D-modules (see [TT], [T3]) to get a holonomic system for the integral of a function 
involving the logarithm of a polynomial. 

2 Annihilators with a parameter 

Let / be a non-constant polynomial in n variables ) with coefficients in 

the field C of the complex numbers. From an algorithmic viewpoint, we assume that the 
coefficients of / belong to a computable field. 

First, we consider formal functions of the form f 8 (logf) k with an indeterminate s. 
More precisely, for a non-negative integer m, we introduce the module 

in 

C(f,m) :=0C[x,f 1 , S ]f(log/) fc 1 

of which f s (\ogf) k are regarded as a free basis over C[x, / _1 ,s]. Then £(f,m) has a 
natural structure of left Z) n [s]-module, which is induced by the action of the derivation 
dj = d/dxj defined by, for a e C[x, s], 

^•{ar(log/) fc } = Q^ + sa/- 1 !^) f'Qoef) k + kar 1 ^-f'Oo S f) k - 1 U = l,...,n) 
if k > 1 and 

«,(„/.)= (g+^-ig)/- (,• = !,...,„). 

In view of this action, it is easy to see that £(/, m)/£(f, m — 1) is isomorphic to C(f, 0) = 
C[x, f~ x , s]f s as a left -D n [s]-module. 
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Now consider the left _D n [s]-submodule 

V(f, m):=D n [s]f s + ..- + D n [s](f s (log ZD 
of C(f,m). Our purpose is to determine the annihilator module 

{m 
P = (P , Pi, ... , P m ) G D n [s] m+1 | (log ft) = 

fc=0 

and the annihilator ideal 

Ann Dn[s] f (log /)™ := {P G D n [s] \ P(f s (\ogf) m ) = 0} . 

Note that there are isomorphisms 

P(/,m) ~ P n [ S ] m+1 /Ann D „ w (/ s , . . . , /'(log/)" 1 ), 
D n [s}(f s (log f) m ) ~ D„[ S ]/Ann DnW (r(log/) m ). 

Now let us regard f s \ogf as a multi- valued analytic function in (x, s) on {(x, s) G 
C™ +1 | /(x) ^ 0}. 

Lemma 1 Let f G C[x] be a non-constant polynomial. Then for di(x) G C[x], 

m 

^(xxiog/y = o 

i=0 

ZioWs as analytic function if and only if (Zj(x, s) = /or a// i. 

Proof: We argue by induction on m. Let xq G C n be a non-singular point of the hyper- 
surface f(x) = 0, i.e, assume 

df 

f(xo) = 0, — — (xq) 7^ for some i with 1 < i < n. 

CXj 

In view of the uniqueness of analytic continuation, we have only to show that each ai(x) 
vanishes near xq. Hence we may suppose that ai(x) are analytic near xq = and f(x) = x\. 
That is, 

a (x) + ai(x) logxi H h a m (x)(logXi) m = (2) 

holds on a neighborhood U of 0. Fix a point x = (x±, . . . , x n ) in U such that x± ^ 0. By 
analytic continuation along a circle (e'^~^ t xi, X2, ■ ■ ■ , x n ) with < t < 2tt, the identity (T5]) 
is transformed to 

a (x) + ai(x)(logXi + 2ttv /Z 1) H h a m (x)(logXi + 27rv /Z l) m = 0. 

By subtraction, we get an identity of the form 

b (x) + bx(x) logxi + • • • 6 m _ 1 (x)(logx 1 ) m " 1 = 

with 

b m -i(x) = 2mn\^la m (x). 

From the induction hypothesis it follows that b$(x) = ■ ■ ■ = b m ^i(x) = 0, which implies 
a m (x) = 0. We are done by induction on m. □ 
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3 Computation of the annihilator 

Now let us describe an algorithm for computing the annihilator of f s (log f) m . 

Algorithm 1 Input: a non-constant polynomial / in the variables x = (x\, . . . , x n ) with 
coefficients in a computable subfield of C, a non-negative integer m. 

(1) Let G = {Px(s), . . . , Pfe(s)} be a generating set of the left ideal Ann Dn [ s ]/ s := 
{P(s) G D n [s] | P(s)f s = 0} by using an algorithm of [9] or [3] (see also [5]). 

(2) Let e = (1,0, ... ,0), ... ,e m = (0, . . . , 0, 1) be the canonical unit vectors of C m+1 . 
For each i = 1, . . . , k and j — 0, 1, . . . , m, set 

Output: G' := {Pi( s ) 0) I 1 < * < fc > < J < m ) generates Ann DnH (/ s , . . . , / s (log/) m ). 

Algorithm 2 Input: a non-constant polynomial / in the variables x = (xi, . . . , x n ) with 
coefficients in a computable subfield of C, a non-negative integer m. 

(1) Let G' be the output of Algorithm [TJ 

(2) Compute a Grobner base G" of the module generated by G' with respect to a term 
order -< for (D n [s]) m+1 such that Mej -< M'ek for any monomial M and M' if k < j. 
Let Go be the set of the last component of each element of G". 

Output: G generates Ann Dn ^f s (log f) m . 

Lemma 2 Let I be a left ideal of D n [s] generated by {Pi(s) , . . . , Pk(s)} . For P(s) G 
-D n [s], and j G N, set 

Then the left D n [s]-submodule of (D n [s]) m+1 which is generated by {P(s)v> \ P(s) G I, < 
j < m} coincides with the one which is generated by {P»(s)^ \ < i < k, < j < m} 
for any integer m > 0. 

Proof: Let J\f be the left _D n [s]-module generated by {Pj(s)^ \ < i < k, < j < m} 
and P(s) be a nonzero element of /. Then there exist 

mi 
1=0 

such that P(s) = X/t=i Qj( s )-^( s )- Then we have 

^) (j) = EE^(^( s )) (J) - 

i=l /=0 
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Hence we have only to show that (s l Pi(s))^ belongs to Af . This can be done as follows: 




where (f) M := 1(1 - 1) •••(/- \i + 1). □ 

Theorem 1 The output of Algorithm^ coincides with AnnD n [ s ](/ s , . . . , / s (log f) m ). 

Proof: Let P(s) belong to AnnD n [ s ]/ s . Differentiating the equation P(s)f s = with 
respect to s, we get 

for < j < m. This shows that each Pi(s)v> annihilates (f s , . . . , f s (\og f) m ). 

Set M := Ann Dn [ s ](f s , . . . , f s (\ogf) m ). Let Af be the left D n [s]-module generated by 
the output G' of Algorithm [TJ The argument above shows that Af is a left .D n -submodule 
of AA. Hence we have only to prove Af — AA. For this purpose let Afj be the left 
_D n [s]-module generated by {Pj(s)^ \ 1 < i < k, < u < j} and set 

Mj := {{Q , Qi, . . . , Q m ) e M \ Q v = if v > j}. 
Let Q(s) = (Qo(s), . . . , Qj(s), 0, . . . , 0) be an element of Aij. Then 

J2Q u (s)(f(hgfy)=0 

holds. In view of the action of D n [s] on£(/,j) noted in Section 1, this implies Qj(s)f s = 0. 
Hence Qj(s)^ belongs to Afj by Lemma [2j It is easy to see that Q(s) - Qj(s) U) belongs 
to AAj~\. This means AAj = Afj-\-M.j-\ for 1 < j < m. Then we can show that Afj = A4j 
holds for 1 < j < m by induction on m noting Afo = Aio. □ 

Remark 1 If / is weighted homogeneous, i.e., if there exist rational numbers Wi such that 
Y^i=i w i x i^i(f) = /) then D n [s]f s (log f) m is isomorphic to C(f,m) as left _D n [s]-module. 
That is, we have an isomorphism 

D n [s}/ Ann Dn[s] f (log f) m * (D n [s]) m+1 / Aim Dn[s] (r , / s (log/) m ). 

of left _D n [s]-module In fact, this follows from the relations 

- s) (f(hgf) k ) = kf s (\ogff- 1 (k > 1). 
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4 Specialization of the parameter 

Let us fix a complex number A. (From an algorithmic view point, we assume A lies in a 
computable subfield of the field C.) We set 

m 

C(f,m,X) :=^C[xJ- 1 ]f x 0ogf) i f 

where f x (\ogf) k are regarded as a free basis over C[x, Substituting A for s gives 

C(f,m, A) a natural structure of left -D n -module. In fact, one has 

^•{a/ A (log f) k } (J£- + W" 1 ^) />g ff + kaf~ l ^f \\og ft' 1 (j = 1, . . . , n) 
for k > 1 and 

%<./») = (£ + *>r'g)/» w = i,....») 

with a G C[x, This implies that £(/, m, \)/C(f, m—1, A) is isomorphic to £(/, 0, A) = 
D n f x as a left D n -module. It follows that C(f, m, A) is holonomic since so is D n f x as was 
proved by Bernstein [1]. 
Set 

V(f, m, A) := D n f x + ■■■ + D n (f x (log f) m ) 
We define the annihilators of . . . , f x (\ogf) m ) and of f x (\ogf) m to be 

m 

Ann D „ (/\ . . . , / A (log f) k ) := {P = (P ,Pi, P m ) G (D n ) m+1 \ ^ P fe (/ A (log f) k ) = 0}, 

fc=0 

Ann Dn / A (log/) m := {P E D n \ P(f x (\ogf) m ) = 0} 

respectively. Then we have isomorphisms 

V(f, m, A) ~ (D n ) m+1 /Ann Dn (f x , . . . , / A (log ft) 
D n {f\\ogf) m ) ~ D n /Aim Dn {f x ()o S f) m ). 

In the sequel, we need information on the integral roots of the Bernstein- Sato poly- 
nomial or the h-function of /, which is, by definition, the monic polynomial bf(s) of the 
least degree such that a formal functional equation 

p( S )r +i = b f (s)f (3) 

holds with some P(s) G -D„[s]. The existence of such a functional equations was proved 
by Bernstein [TJ. It was proved by Kashiwara [4J that the roots of 6/(s) = are negative 
rational numbers. An algorithm to compute 6/(s) and an associated operator P(s) was 
given in [8]. The following proposition generalizes a result of Kashiwara [U Proposition 
6.2]: 
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Theorem 2 Let bf(s) be the Bernstein- Sato polynomial of f , i.e., a polynomial in s of 
the least degree such that P(s)f s+1 = bf(s)f s holds with a P(s) G .D n [s]. Let X be a 
complex number such that bf(X — v) ^ for any positive integer v . Then we have 

Ann Dn (/\...,/ A (log/) fc ) = {P(A) | P(s) G Ann DnH (r,...,r(log/) fc ), 
An nA J A (log/r = {PW I P(s) G Ann DnH r(log/r}. 

Proof: We have only to show the first equality. Assume that ^^L Pfc(/ A (log f) k ) = 
holds with Pk G D n . Then there exist non- negative integer I > and polynomials 
a k (x, s) G C[x, s] such that 

m m 

Pk(f (log f) k ) = (s-X)Y,a k (x,s)f s - l (logf) k . 

k=0 k=0 

By using the functional equation d3J), we can find an operator Q(s) G D n [s] such that 

b f (s-l)---b f (s-l)f- l = Q(s)r. 

In view of the action of D n [s] on C(f, to), there exist a' k (x, s) G C[x, s] and a non-negative 
integers Zi such that 

m— 1 

6/ ( a _ i) . . . 6/ ( a _ or-<(io g /r = g( a ){/-(iog /n + s )f s ' h ( lo g ff- 

fc=0 

Proceeding inductively, we conclude that there exist a polynomial b(s) G C[s] which is a 
product (possibly with multiplicities) of bf(s-j) with j > 1 and operators Qk(s) G D n [s] 
such that 

m m 

6(a) £ s)r- z (log /) fc = £ Q*(a){/'(log /) fc }- 

fc=0 fc=0 

Hence 

m m 

P(s) := b(s) Y P^k ~(s-X)J2 Qk{s)e k 

k=0 k=0 

belongs to Ann DnH (f, . . . , / s (log f) k ) and 6(A) YJk=o p k e k = -P(A). This completes the 
proof since 6(A) ^ by the assumption. □ 

If 6/(A — v) — for some positive integer i/, then set vq := max{i/ G Z | 6/(A — u) = 0} 
and Ao := A — Uq. Then Ao satisfies the condition of Theorem [2j Then for (P , . . . , P m ) G 
(P„) m+1 , we have 

(P ,...,P m )GAnn Dn (/ A ,...,/ A (log/) fe ) 

& (PoF°, P m D G Ann Dn (f Xo , • • • , / Ao (log /)*) 
4» (Po, • • • , P m ) G Ann D „(/ A °, . . . , / A °(log/) fc ) : f°. 

The module quotient in the last line can be obtained by computing the module intersection 
or else by syzygy computation. Now let us describe two algorithms for module quotient in 



7 



general. First, let us define the componentwise product of two elements P = (P , . . . , P m ) 
and Q = (Q , • • • , Qm) of {D n ) m+1 to be PQ := (P Q , P m Q rn ). Let N be a left 
-D n -submodule of (D n ) m+1 and P be a nonzero element of (D n ) m+1 . Then the module 
quotient N : P is defined to be 

N:P:={Qe (D n ) m+1 | QP G iV}, 

which is a left _D n -sub module of (D n ) m+1 . 

Algorithm 3 Input: A set G\ of generators of a left -D n -submodule N of (D n ) m+1 and 
a non-zero element P = (P , Pi, • • • , Pm) of (Ai) m+1 - 

(1) Introducing a new variable t, compute a Grobner base G 2 of the left _D n [t]-module 
of (D n [t]) m+1 which is generated by {(1 - t)P k e k | < k < m} U {tQ \ Q G d} 
with respect to a term order -< such that x a d /3 ej -< te k for any j, k G {0, 1, . . . , m} 
and a, (3 E N n . 

(2) G 3 := G 2 n (,D n ) m+1 . 

(3) G 4 := {Q/P | Q G G 3 }, where Q/P denotes the element in {D n ) m+1 such that 
(Q/P)P = Q in the sense of compnentwise product. 

Output: G 4 generates the module quotient iV : P. 

In fact, we can show in the same way as in the commutative case that G3 generates the 
left module Nn(D n ) m+1 P. In particular, for each Q G G 3 , there exists Q' G (D n ) m+1 such 
that Q = Q'P. Let us denote this Q' by Q/P. Then Q' belongs to the quotient module 
N : P. Conversely, if Q' belongs to : P, then Q'P belongs to N n (D n ) m+l P. Hence 
Q' belongs to the module generated by G 4 . The correctness of the following algorithm 
should be clear: 

Algorithm 4 Input: A set G\ = {Qi, . . . , Q k } of generators of a left -D n -submodule 
of (D n ) m+1 and a non-zero element P = (P , . . . , P m ) of (P>„) m+1 . 

(1) Compute a set G 2 of generators of the syzygy module 

m m+k 

S := {(So, S±, . . . , S m , S m+ i, . . . , S m+k ) G (D Tt ) m+fc+1 I SjPjCj+ S m+ jQj = 0} 

j=0 j=m+l 

via a Groebner base of the module generated by Pjej (0 < j < m) and Qj (1 < j < 
k). 

(2) Let G 3 be the set of the first m+1 components of the elements of G 2 - 
Output: G3 generates the module quotient : P. 

Summed up, the annihilators for (/ A ((log/) fc ) <fc< m and (/ A (log/) m are computed as 
follows: 
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Algorithm 5 Input: a non-constant polynomial / in the variables x = (x\, . . . , x n ) with 
coefficients in a computable subfield of C, a number A which belongs to a computable 
subfield of C, a non-negative integer m. 

(1) Compute a set Gi of generators of Ania Dn ^(f s , . . . , / s ((log/) fe ) by Algorithm [TJ 

(2) Compute the (global) Bernstein-Sato polynomial 6/(s) of / by using one of the 
algorithms in [8], [9], [3] or their modifications. 

(3) Let Uq be the largest positive integer v such that bf (A — i>) = if there are any such 
v. If there are no positive integer v such that 6/(A — i/) = 0, then set 2/0 = 0. 

(4) Set Ao := A — uq and G 2 '■= Gi|s=a (substitute Ao for s in each element of G\). 

(5) If uq > 0, then let G3 be a set of generators of the module quotient (G 2 ) '■ f u ° = 
{G 2 ) '■ (f v °, ■ ■ ■ , f u °), where (G 2 ) denotes the left module generated by G 2 - 

(6) If i/ = 0, then set G 3 := G 2 . 

(7) Compute a Grobner base G4 of the module generated by G3 with respect to a term 
order -< for (D n ) m+1 such that Mej -< M'ej- for any monomial M and M' if k < j. 
Let G5 be the set of the last component of each element of G4. 

Output: G 3 generates Arni Dn (f x , f x (\ogf) m ); G 5 generates Ann Dn / A (log f) m . 

Remark 2 In step (3) of the algorithm above, we need only integer roots of the 6-function. 
Hence one can employ a method described in [6] to determine all the integer roots of the 
fe-function efficiently without computing the whole 6-function. 

5 Implementation and examples 

We have implemented the algorithms in a computer algebra system Risa/Asir [7], which 
is capable of Groebner base computation of modules over the ring of differential operators 
as well as over the ring of polynomials. 

Example 1 (one dimensional case) Let / be a square-free polynomial in one variable 
x with complex coefficients. Since AnnDj[ s ]/ s is generated by fd — sf, the annihilator 
module Ann Dl [ s ](/ s , . . . , / s (log/) m ) is generated by m + 1 elements 

(fd x - sf, 0, ■ ■ • , 0), (-/', fd x - sf, 0, • ■ ■ , 0), • ■ ■ , (0, ■ • • , 0, -mf, fd x - sf) 

with d x = d/dx and /' = d(f). 

Since the Bernstein-Sato polynomial of / is s + 1, the substitution s = A gives genera- 
tors of Ann Dl (/ A , . . . , / A (log f) m ) if A ^ 0, 1, 2, .... In particular, Ann Dl (f~\ f-\hg f) m ) 
is generated by 

(3 x /,0,...,0), (-/',0*/,O,...,O), (0,... ,0,-mf,dJ). 
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In view of Algorithm^ we can verify that Ann^l, (log f) m ) = Ann^/ l ,...,f 1 ()og 
f is generated by 

(d x , 0, • ■ ■ , 0), (-/', fd t , 0, • • • , 0), ■ ■ • , (0, ■ • • , 0, -mf, fd x ). 

Explicit generators of Ann^^log f) m for m > 1 would be complicated: For example, if 
f = x 3 — x and m — 1, Algorithm [5] gives generators 

(3x 5 - 4a; 3 + x)<9 2 + (3x 4 + 1)9 X , 

(a; 3 - x)<9 3 + (-3a; 4 + 9a; 2 - 2)d 2 x + (-3a; 3 + 3x)d x 

of Ann^j log/, which is not generated by a single element. 

Example 2 Set / = x 2 y 2 + z 2 with n = 2 and (x%, x%, X3) = (x, y, z), d x — d/dx and so 
on. First Aian D3 ^f s is generated by 

-xd x + yd y , yd y + zd z - 2s, 
zd x - y 2 xd z , zd y - yx 2 d z , 
-zd 2 x + y 3 d z d y + y 2 d z . 

Since the Bernstein-Sato polynomial of / is bf(s) — (s + l) 3 (2s + 3), the substitution 
s = — 1 gives a set of generators of Arna^f' 1 log/. Then by ideal quotient computation 
we get a set of generators 

-xd x + ydy, -zd x + y 2 xd z , 
d 2 + x 2 dl d 2 + y 2 d 2 z , 
-zdy + yx 2 d z , dydl - zyd 3 , 
-d 2 y d x + zxd 3 z , yd y d x + zd z d x 
yd z d y + zd 2 z + d z , 
yd 2 + zd z d y + d y , 
zd y d x + zyxd 2 - yxd z , 
d 2 d 2 x + z 2 &\ + 2zd 3 

of Anno 3 log/. Let us consider the integral 

u(t) := ! e -t(* 2 +y 2 +z 2 ) \ og (x 2 y 2 + z 2 ) dxdydz, 

which is well-defined for t > 0. Then u(t) satisfies ordinary differential equations 

P lU (t) = P 2 u{t) = 

with 

P l = t 3 8f + (2t 4 + I7t 2 )d? + (32t 3 + 80t)<9 3 

+ (-At 4 + lUt 2 + 100)<9 2 + (-28t 3 + 192t)d t - 36t 2 + 48, 

p 2 = t 3 8? + (3t 4 + Ut 2 )d 3 + (2t 5 + 35£ 3 + 52t)d 2 
+ (Ut A + 102t 2 + 48)<9 4 + 18t 3 + 66t. 
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6 (Appendix) Holonomic systems 



Let us present a precise definition of holonomicity. We define the total or the (1, l)-order 
of nonzero P E D n to be 

ord(i ] i ) (P) := max{|a| + = «i H h a n + (3\ H h /3 n | a a ^ ^ 0}. 

We set ord„,(0) := — oo. This induces the filtration 

F k {D n ) := {P E D n | ord (lil) (P) < k} [k E Z) 

on the ring D n . Let M be a left _D n -module and {F k (M)} k£ z be a good (1, l)-filtration. 
This means the following properties: 

(1) every F k (M) is a finite dimensional vector space over C; 

(2) F fc (M) C F fc+ i(M) for all fceZ; 

(3) |jF fe (M) = M; 

fcGZ 

(4) F t (D n )F k (M) C F i+fe (M) for all i,k E Z; 

(5) there exists fc^Z such that F k (M) = for k <k\] 

(6) there exists k 2 E Z such that Fi(D n )F k (M) = F i+k (M) for > fc 2 . 

Then there exists a polynomial in /c such that dimc-Ffc(M) = p(k) for sufficiently large 
A;. The degree of p(k) does not depend on the choice of a good (1, l)-filtration of M and 
is called the dimension of the module M, which we denote by d(M). It was proved by 
Bernstein pQ that d(M) > n if M ^ 0. The following definition is due to Bernstein pQ: 

Definition 1 A finitely generated left _D n -module M is called a holonomic system if 
d(M) < n. We also call a left ideal I of D n to be a holonomic ideal, by abuse of 
terminology, if the left _D n -module D n /I is holonomic. 

Note that <i(M) < n is equivalent to d(M) = n or M = in view of the Bernstein 
inequality stated above. The dimension d(M) can be computed as the degree of the 
Hilbert function from a Grobner base with respect to a term order which is compatible 
with the total degree. 

Holonomicity is preserved by operations such as sum, product, restriction to affine 
subvarieties, and integration with respect to some of the variables (cf. [I], [2]) and they 
are computable (see e.g., [10]). Let R n := C(x){d\, . . . ,d n ) be the ring of differential 
operators with rational function coefficients. A D n -module M is said to be of finite rank 
and the dimension is called the rank of M, if R n M is a finite dimensional vector space over 
C(x). A holonomic -D„-module M is of finite rank but the converse is not true in general. 
Note that there is an algorithm for a given D„-module M of finite rankjo construct a 
holonomic D n -module M and a surjective D n -homomorphism of M to M (|12j.[T4"j). If 
M is a system of differential equations of finite rank for an analytic function u, then we 
have an isomorphism ~D n / Aim.£, n u. 
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Example 3 Set / = x 2 y 2 + z 2 and consider the function / l . It is easy to see that the 
operators 

f9y + % = {x2yl + ^ + 2x2yi 

fdz + % = ^ + ^ + 2Z 

annihilate f~ l . Let J be the left ideal generated by these three operators, a 'naive' 
annihilator. Then the Hilbert function of the D 3 /J is 

1 , 5 1 , 4 7 3 5 2 43 , 

which means that the degree of D 3 / J is 5 and hence D 3 /J is not holonomic although it 
is of rank one. The true annihilator / of f~ l is generated by 

3z 2 d 2 x - 2y 3 d z d y - 2y 2 d z , 
3z 2 d y — 2yx 2 d z , 
3z 2 d x — 2y 2 xd z , 
Sydy + 2zd z + 6, 
-xd x + yd y 

and the Hilbert function of D 3 /I is 

7 ,3 3 ,2 4 3, , 

which implies that D3/I is holonomic. The Hilbert function of D 3 / Am\£, 3 log / is 

2A; 3 + -k 2 + -k - 1. 
2 2 
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